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Abstract 

We study the residual gauge freedom within the null quasi-spherical (NQS) 
gauge for spacetimes admitting an expanding shear-free null foliation. By 
constructing the most general NQS coordinates subordinate to such a folia- 
tion, we obtain both a clear picture of the geometric nature of the residual 
coordinate freedom, and an explicit construction of nontrivial NQS metrics 
representing some well-known spacetimes, such as Schwarzschild, accelerated 
Minkowski, and Robinson- Trautman. These examples will be useful in test- 
ing numerical evolution codes. The geometric gauge freedom consists of an 
arbitrary boost and rotation at each coordinate sphere — and this freedom 
may be used to normalise the coordinate to an "inertial" frame. 

04.20,04.30 



Typeset using REVTeX 



* from March 1997: robertb@ise.canberra.edu.au 



1 



I. INTRODUCTION 



The recently introduced [|T]-^ null quasi-spherical (NQS) coordinate condition provides 
a new approach to the study of the Einstein equations in exterior regions admitting an 
expanding null foliation. The NQS gauge is described by the metric ansatz 

ds%Qs = -2u dz{dr + v dz) + [rd'd + jS^dr + 7^^)^ + (r sin M^p + 13'^dr + 7^dz)^ (1) 

where {'d, ip) are the usual polar coordinates on S'^, m > and v are real-valued functions, 
and 

(3 = (3^dii + (3''csc'dd^, j = j'd^ + j\sc^d^, (2) 

may be considered either as vectors tangent to the spheres {z, r) = const or, using a complex 
formalism, as spin-1 fields. 

The advantages of the gauge, and its generality, are discussed in [Q. The purpose of 
this paper is to analyse the gauge freedoms remaining within the NQS gauge condition, for 
the class of spacetimes admitting a shear-free (/3 = 0) null foliation. We explicitly describe 
the construction of such foliations in Schwarzschild, Minkowski and Robinson- Trautman 
spacetimes. 

The examples will also be useful as test data for numerical solvers, since they involve 
arbitrary functions but are still simple to describe explicitly. This remark applies both 
to characteristic and 3+1 codes — the class of boosted Schwarzschild metrics should be 
particularly appropriate as test data. 

The NQS gauge is best understood by comparison with other popular conditions used 
to describe the metric on a null foliation, due to Bondi |^ and Newman and Unti 0. The 
Newman-Unti radial coordinate r is determined by a choice of affine parameter along each 
of the null generators; the Bondi radius is defined by the condition that the spatial volume 
form sm{}d'd A dip have length (r^ sin'i?)"^. In both cases the angular coordinates (i?, (p) are 
transported along the null generators. Both coordinate systems are determined by labelling 
and normalisation conditions at just one transverse S"^ in a null hypersurface and therefore 
have gauge freedom corresponding to functions on a single S*^ (in each null hypersurface). 

By contrast, the NQS radial function r has level sets isometric to standard spheres of 
radius r. Although it is possible to then determine the angular coordinates {■&, ip) as labelling 
the outgoing null generators (as in p, 0), it seems more geometrically natural to use the 
("i?, (/?) determined by the isometry with S"^ with metric r'^^d'd'^ + sit? "d dtpP'). 

Since the metric spheres at each radius are not unique within the null hypersurface 
(at least, this is the case for the standard Minkowski null cone), there is an additional 
coordinate freedom within the NQS gauge, consisting of a choice of Lorentz transformation 
at each sphere. This freedom does not have an analogue in the Bondi and Newman-Unti 
gauges. 

The vector field (3 is referred to as the shear, that this terminology does not conflict 
with the accepted usage of "shear" is seen by noting that the (usual) shear of the null 
generator i = dr — r~^j3 of the NQS metric Eqn. (0) is given in the Newman-Penrose 
notation by aNP = t^^S/?, where S is the eth operator on the standard S'^ and we identify 
P ~ -^{P^ — iP'^) with a spin-1 field on S"^. Consequently, vanishing shear vector /? implies 
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vanishing a^p] conversely if a^p = then (3 consists purely of £ = 1 spin-1 spherical 
harmonics The role played by the the ^ = 1 spin-1 spherical harmonics is discussed 
in greater detail in the following section. In the Appendix we show that any shear-free, 
expanding and twist-free metric admits NQS coordinates with /3 = — this was shown in 
||Tl| , p^ for vacuum metrics. 

The metric form (0) with /? = 0, when restricted to a coordinate null hypersurface C, 
becomes 

dsl = r\d^^ + sm'^^dip'^). 

By identifying C with the future null cone at the origin in M'^'^, we can see that this form 
is invariant under the Lorentz group 5*00(3, 1) — and the Lorentz transformation may also 
vary with r, since invariance only requires that each quasi-sphere r = const, is mapped 
isometrically. Thus, our main idea is to use explicit representations of the Lorentz group 
acting on the standard null cone Co = {t = \x\} in Minkowski space M^'^ to describe the 
general transformation leaving the form ds^. invariant. 

Note that the problem of finding general quasi-spherical foliations of a null hypersurface 
which is not shear-free and expanding is considerably more difficult, since the explicit model 
of the standard cone and its associated Lorentz deformations is no longer available. However, 
linearisation arguments suggest strongly that the gauge freedoms of the shear-free case are 
mirrored in the more general setting, provided the shear is not too large. Thus we expect 
that the description here of the shear-free NQS freedom will provide some insight into the 
more general case. 

At least in the shear-free case, we will show that the NQS condition has gauge freedom 
consisting of an SOo{S, l)-valued function of the radius (on each null hypersurface); this is 
functionally less rigid than the Bondi and NU gauges, since it has freedom in r which is 
lacking in these gauges. This Lorentz transformation freedom may be viewed as providing 
a choice of "inertial frame" normalisation at each radius, and may be used to normalise 
certain of the remaining metric coefficients, as described below. 

This interpretation is supported by a comparison ^ between the Robinson- Trautman 
metrics and the NU form of the Minkowski metric in coordinates using null cones 
with base point describing a timelike curve. This comparison may also be made in the 
NQS coordinates, and supports both the interpretation of Robinson- Trautman spacetimes 
as describing an accelerated black hole rapidly settling down to a Schwarzschild black hole 
in uniform motion, and the interpretation of the NQS freedom as representing a choice of 
reference frame at each radius and time. 

In section 2 we study the metric and NQS freedom of the model cone Co, by constructing 
the most general quasi-spherical (QS) foliation of Cq. The resulting metric has shear vector f3 
consisting solely of £ = 1 spherical harmonics, and we show conversely that any null surface 
with such shear vector is gauge-equivalent to the standard cone. In section 3 we describe 
the metric in general NQS coordinates on a spacetime admitting a shear-free null foliation. 
Section 4 describes the application of these results to the specific examples of Schwarzschild, 
Minkowski and Robinson- Trautman spacetimes. Basic results on shear-free expanding null 
hypersurfaces are collected in the Appendix. 

The computations are presented in slightly more detail than is strictly necessary, in order 
to facilitate the use of the example NQS metrics in benchmarking numerical codes, and in 
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the interpretation of general NQS numerical results. 

II. MODEL CONE Cq 

Let Co = {{x,t) G M:''^ : t = \x\} be the standard future null cone based at the origin in 
Minkowski space. We may use x G M'^ as a coordinate on Co; instead a polar representation 

x — r9, r — \x\, 9 — x/r — (xi/r) & S'^ (3) 

will be very useful, where we identify = {x G : |a;| = 1} and we use the direction 
cosines 9 — (^i, ^2, ^3), \9\ — 1 to parameterise 5"^. The polar coordinates on Co will usually 
be denoted by (r, 9) or (p, (). 

The parameterisation (^j) of 5^ leads to a representation of tangent vector fields to S"^ 
as 3- vector fields on the unit sphere 5*^ C which are tangent to S^. This will prove more 
convenient than using the polar coordinate basis {d^, d^). Thus, a vector field Y — Y{9; A), 
depending on ^ G 5"^ and other parameters A (eg A = {z,r)), may be represented as the 
3- vector Y = (Yi) satisfying e^Y{9, A) = 9iYi = 0. 

Throughout we use latin indices i, j, ■ ■ ■ with range 1, . . . , 3 and the summation convention 
on repeated indices, not necessarily raised and lowered. 

The Minkowski metric induces the rank-2 degenerate bilinear form 

dsl^r'\d9\' = r'j2id9,r (4) 

i=l 

on Cq. Note that \d9\'^ is the standard metric on S'^, and \d9\'^ = Yl,ij=i ^ijdxidxj, where 6 
is the projection matrix 

e = / - 99'^, Qij = % - 9i9j, (5) 

and 9^ represents the transpose (row) vector. 

A quasi-sphere of radius r G M.'^ in Co is an orientation preserving embedding '■ ^ Co 
such that 

Kidsl)=r'\d9\', 

and we say $r is a quasi- spherical map. 

Every quasi-sphere in Cq is determined by a unique time and space orientation preserving 
Lorentz transformation L G 5*00(3, 1) via the composition $r = L o with the inclusion 
V : 52 ^ Co, 9 ^ [r9^, rf, 

e Co, (6) 

and we now exploit this basic description. 

The Lorentz transformation L G 500(3, 1) admits a unique boost/rotation decomposition 
L = nB with 



^r:9eS' 



r9 
r 



gCo k 



r9 
r 
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n = n{R) = 

B = B{w) = 



R 
1 

W w 
b 



R e S0{3), 



w e 



(7) 
(8) 



where b — y^l + \w\'^ > 1 and W is the 3x3 matrix 

1 



b+1 



WW 



(9) 



It is easily checked that B preserves the Minkowski metric rj, r] — B^rjB. The quasi-spherical 
map $r iiiay be described explicitly by 



MO) 



' r-e' 




r 





rR9 + r(l + ) Rw 
r(b + e'^w) 



(10) 



where 9^w = w^9 is the usual inner product between column 3-vcctors. In terms of the 
rectangular W^'^ parameterisation of the target Cq we have ^r{9) = [p{r,9)({9)^ , p{r,9)]'^ , 
where we define 



f{9) = f{9;w) 
p(r,9) = p{r,9;w) 

a9)^ao;w,R) 



= b + 9'^w^ Vl + + 9^w, 
= rf{9) 

= r'R + (i + 6f) ^) = r'Ri^ + WO)- 



(11) 

(12) 
(13) 



Here and elsewhere we adopt the convention that (p, () denotes polar coordinates on the 
range (target) Cq of which leads to the description $^(6') = {p{r,9)X{0)) of in polar 
coordinates. The metric ds'^^ on the target cone Cq in polar coordinates (p, Q is just p'^\dQ^ 
and we may verify the quasi-spherical condition ^*{ds^^) — r^\d9\'^ by direct computation 
as follows. 

We define the angular gradient operator Dq = {Dq.) as the projection tangent to the 
unit sphere 5"^ of the ordinary gradient in W?. Exphcitly, let h(.{x) := h{x/\x\), x e 
be the homogeneous degree extension of any h e C^{S'^) and define De^h — dhe/dxi. It 
follows that 9iDg^h = by the homogeneity condition, and Do^9j — Qij. Then ^*{p'^\d(\'^) — 
p{r,9)%ed9\'^, where {C,e)ij = Dq.Q with 



Dg.Ci = Dg^ {{b + 9^w)-'Rik {Wk + Wkl9l)) 

= r^Rik {-r^{Wk + Wkl9i)QjrnWm + W^fcm©mj) 
— f Rik-^km^ mj ■) 

and we have introduced the very useful matrix 

A := 7 - /-^ {9 + {b + l)~^w) w^. 



(14) 



(15) 



By exercising a certain amount of care we may convert to a matrix notation. We adopt 
the convention that 3- vector quantities such as 9i, d9i, Wi, are to be treated as column vectors 
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(of 1-forms, functions, etc.) and that row vectors usually will be indicated by the transpose 
notation, except that we regard Dg, Z)^ and their associated gradients as row vectors. 
For example, with these conventions the computation ( [I^ ) may be summarised as 



= f-'RAQ 

and the chain rule for h G C°°(S'^, M) appears as 

dh = h g dO 
d{hoQ = h^^C,ede 
{hoC)^s = r^h^RAQ, 

where both sides of the final identity are row vectors. 
The matrix A satisfies several useful identities, 



(16) 



A^A 
GA^AQ 

A-' 



I - f-\we^ + ew^) 
e, 

I+{9 + {b+V-' 



w )w 



(17) 
(18) 
(19) 



which may be used to show that $r is quasi-spherical: 



r'^P\r^RAdd\^ = r^de^A^R^RAde 
= r^de^A^Ade = r^\de\'^, 

since R^R = I and 9^d9 = 0, so d9 = Qd9. This confirms that $r is quasi-spherical and 
also conformal, since ^*{\d(\'^) = f~'^\d9\'^. Note that from Eqs. (|T3D ,(|T9D we have 



C(^) = r'RA-'''9, 



(20) 



where A~^'^ is the inverse transpose matrix. 

Having described a single quasi-sphere, we may now consider the effects of r-dependence: 
a map $ : Cq — Cq is said to be quasi-spherical if $r = ° V is a quasi-spherical map for 
each r > 0, and if r i— $r is at least continuously different iable in r — although for 
simplicity we shall consider only smooth maps. Equivalently, $(r, ^) = L{r) o ir{9), where 
the Lorentz transformations L{r) = TZ{R{r))B{w{r)) are described by boost and rotation 
maps w e C°°(R+,M3), R G C°°(M+, ^0(3)). Note that we do not require that $ be a 
diffeomorphism, and this does not follow from the condition ^*{ds'^^) = \d9\'^ since ds"^^ is 
degenerate. Note also that when defined, is not usually quasi-spherical; neither is the 
composition of two quasi-spherical maps usually again quasi-spherical. 

The general quasi-spherical map $ : Co — Co is thus described by 



^r,9) = n{R{r))B{w{r)) 



where the functions p(r, 6'), C('", 6') are defined by Eqs. (^),(|T^) with w,R now depending 
on r. The pull-back metric is then 



'r9' 




p{r,9)Cir 


r 




p{r,9) 
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<^*{dsl)=p{r,9)Xed9 + Crdr\' 

= p%ed9\^ + 2p'de'^CKrdr + p'lCrprfr^, 

where C,r denotes the column vector of partial derivatives ^C,. Note that we are taking the 
liberty of using p, ^ to denote both the coordinates p, C on Cq and their pullbacks p(r, 6) = 
$*(p), ({r, 6) = $*(C); which are functions determining the map $ — this ambiguity should 
not lead to any serious confusion. 

We have already checked that p'^\(^gd6\'^ = r'^\d6\'^, and we compute from Eqn. 



fA-'R-'C,r = A-^R-^R^A-^'^e + A-^-^{A'^'^)e - f-^fM^Ay^O. (21) 

The first term on the right hand side is simplified by introducing the antisymmetric matrix 

■= WR-^R,rW, (22) 
and equals (bearing in mind the formulas = I+b^^6w'^ , W~~^ = I—b^^{b+l)^^ww'^) 

After some computation, we find that the second and third terms of Eqn. (^) may be 
combined into 

'6+1 ^ ' ' 

where a x b = {eijkajb^) is the usual cross product in R'^. Defining the 3-vectors Si,ti 
(depending on r but independent of 9) 

Si := *Si + j;^w^r X w (23) 
h ■■= b-^Siw + W-^w,r (24) 

where *S := (^eijkSjk), and the (r, 6')-dependent vector P 

P .- rBti + re X si, (25) 

we have the identity 

pC,r = RA(3. (26) 
This may be used to simplify ^*{dsQ^^): 

p'dO^CfiCrdr = rde^ A^R^RA(3dr 

= rde'^A^Apdr = rdO'^^dr, 

since 9^ (3 = 0. Similarly we find 

p2|^^^|2 = (3TaTrTrA(3 = \f3\^ 

and it follows that ^*{ds'^^) = \rd6 + (3dr\'^, which is in quasi-spherical form with shear 
vector p. 

To summarise: 
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Proposition 1 Suppose ^ : Co ^ Co is a C°° quasi- spherical map. Then $ satis- 
fies Eqn. ^idj ) for some Lorentz boost parameter w G C°°(M"'",M^) and spatial rotation 
R G C°°(M^, S'0(3)). In the rectangular-polar coordinates {r,6) on Co we have 

<l>*{dsl^) = \rde + (3dr\^, (27) 

where the shear vector [3 is defined in terms of w and R by Eqs. (^), ^^), ^^)- 

The angular vector field (3 consists solely of spin-1 i = 1 spherical harmonics. This follows 
from Eqn. (pSj), since for any constant vector t G M^, the angular vector fields Qt = t — {6^t) 9 
and ^ X t satisfy 

et = grad52(^'^t), (28) 
^ X t = Jgrad52(^^t), (29) 

where the complex structure J : TS"^ TS"^ is defined by anticlockwise rotation with 
respect to the outer normal 6' to S*^ C R'^. Now if ei, 62 is an oriented orthonormal frame on 
S"^ (so Ci X 62 = 9) , we define the spin-1 projection of a vector field X = X^ci + X^e2 by 

X ~ ^=j.(x'-iX') (30) 
(note JX = 9 X X = — X^ei + X^e2 ~ — iO^ ^.nd the operator eth by 

g=^(Ve,-iVeJ, (31) 

where V is the standard covariant derivative on S"^. Note that as so defined, ^ and 9 
are frame dependent, hence the use of ~. Defining the basis vector e = :^(ei —162), we 

could instead write the equality X = + and then consider ^ as the coefficient of the 
representation of X as a section of a spin-1 complex line bundle, with respect to the basis 
vector e. 

l^(j),^pe C\S^,R) then 

5(0 + i7A) = ^(<^,l+ 7^,2 +i(^,i-</',2)), 

where the subscripts (■),„ for a = 1, 2 denote directional derivatives with respect to the basis 
vectors €1,62, and the vector field correspondence is 

grad5.2 — Jgrad52'?/' ~ 5(0 + 1'?/'). 

In particular, for any constant s, t G M'^ we have the correspondence 

et + 9xs ~ g(^^(t-is)), (32) 

and the identity 

As2{9^t) = -2 9'^t (33) 

completes the identification of Bt + x s as a spin-1 i = 1 spherical harmonic. 
We also derive from Eqn. (^) that 

As.{{9^ty-im = -6{{9^tr-l\t\') (34) 



and thus 9^t 9'^s — ^t^s is an £ = 2 spherical harmonic for any s, t G 



P3 



The relations (^3]) , (|2^) between Si,ti and Si,w^r may be inverted, since by direct com- 
putation we have the following lemma. 



8 



Lemma 2 Suppose s,t,w,w,a G and let b := a/1 + anc? W ^ = I — ■^pp^^ww'^ . 
The equations 

t = b~^w X a + W~'^w 

1 ~ v-jS) 

s = a + X "U^ 

are equivalent to 

w = bt + s X w 



a = bW-^s + X t. 



(36) 



Applying this lemma with s = Si, t = ti, w = w^r, a = *Si and assuming (H),(@) gives 
the ordinary differential equations 

w^r = bti + Si X w (37) 
Rr = RW-\b*{W-^si) + ^*(«; X ti))W-\ (38) 

where for any vector t we define the matrix *t = {eijktk)- (So the two star operations 
interchange vectors with antisymmetric matrices.) Consequently we have the following re- 
construction results. 

Proposition 3 Suppose si,ti G C°°(M^,]R'^) are given functions. Given initial conditions 

w{ro) = Wo G , . 

i?(ro) = Ro e 50(3), ^"^^^ 



there is a unique quasi- spherical map $ : Cq — > Cq with parameters w G C°°(]R^, 
R G C°°(M"'", 5*0(3)) satisfying w(ro) = Wq, -R(ro) = Rq and such that the shear vector (5 
satisfies 

j3{r,9) = eti{r) + ex si{r). 



Proof: With given functions of r and b = a/1 + |tf P, Eqn. (|37|) gives an ordinary 

differential equation for w, with initial condition w(ro) = Wq. Since for r G [vq, ri] 

\b{r)ti{r) + si(r) x w{r)\ < 2 (|w(r)| + 1) sup (|ti| + |si|), 

[ro,ri] 

by Gronwall's inequality the solution of Eqn. ( p7| ) is locally bounded and may be continued 
to all r G M"*". Substituting w(r) into Eqn. (|38|) gives an ode for i?(r) G 50(3) with initial 
condition R^tq) = Rq, which similarly has a global solution R G 0°°(M"'", 50(3)). The 
quasi-spherical map defined by the solutions w{r), R{r) via Eqn. ([T0|) has the required shear 
vector P by previous computations. $ is the unique map satisfying the initial conditions at 
tq since any quasi-spherical map may be put into the form (p!OD and the parameters are then 
uniquely determined by the initial value problem (|37| ) ,(|38| ) ,(|39|). QED. 
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Corollary 4 Suppose \rd9 + j3dr\ is a quasi- spherical form on Cq, with shear vector f3 
consisting solely of i = 1 spherical harmonics (ie. f3 may be expressed in the form Eqn. ^Bj)}. 
Then there is a quasi- spherical map $ : Cq ^ Cq, ^{r,6) = (p, C)? such that 

<^*{p^\dC\^) = \rde + f3dr\^ (40) 

Furthermore, $ is unique up to a rigid Lorentz transformation ofCo: if ^ : Cq ^ Cq is any 
map satisfying Eqn. ^4^), then there is Lq G 5*00(3, 1) such that $ = Lq o $. 

Proof: Since f3 is pure i = 1, the i = 1 spherical harmonic coefficient functions Si{r),ti{r) 
are uniquely determined by Eqn. (pSj), and an appropriate quasi-spherical map $ may be 
constructed using Proposition || and initial conditions w(ro) = 0, -R(ro) = / at some radius 
rQ. If $ : Co — > Co also satisfies Eqn. ( |iO|) then $ is quasi-spherical and hence may be 
parameterised by Lorentz transformations L{r), with parameter functions w{r) and R{r). 
Let Wq = w{ro) and -Ro = R{i^o) and let Lq be the corresponding Lorentz transformation. 
Because Ll{r^\de\^) = r^\de\^, the map Lo o $ is also quasi-spherical satisfying Eqn. (^0|), 
and has parameters w{r), R{r) satisfying the initial conditions w(ro) = wq, R^tq) = Rq. 
Since the parameters Si,ti are determined uniquely from /? in Eqn. (^OD, uniqueness of the 
solution of the initial value problem Eqs. (P7D,(PBD,(|55D implies $ = Lq o $ as required. 
QED. 

Note that w{r), R{r) can be computed in terms of w{r), R{r) and Wq,Rq from the 
identity 

L = n{R)B{w) = n{Ro)B{wo)n{R)B{w). 

However it is not true in general that the composition of quasi-spherical maps is again 
quasi-spherical — L in this identity defines a quasi-spherical map only when wq, Rq are 
constant. 



III. DEFORMATION OF SPACETIME METRICS 

We consider now those spacetimes whose metric can be placed in the form 

dslp = -2U dz{dr + Vdz) + \rde + T dz\^, (41) 

where F = T{z,r.,9) is an angular vector field, so F satisfies ^^"^F = 0. We may verify that 
the null congruence defined by the coordinate tangent vector is expanding, shear-free 
and twist-free. This class includes the Schwarzschild, Robinson- Trautman and accelerated 
Minkowski spacetimes, and will be further discussed in the following section. For the present, 
we use the techniques of the previous section to study the effect of quasi-spherical Lorentz 
deformations of metrics of the form (|4TD . 

We regard ds'^gp as defined on (a subset of) M x Co, and then the metric induces the 
standard form r'^\d9\'^ on Co. The form (^) is in fact the most general metric form compatible 
with this property and such that the coordinate z is null (characteristic). Extending previous 
definitions, for any domain C M x ]R+ with coordinates (-z,t), we say that $ : i7 x S*^ — >■ 
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M X Co is quasi- spherical if the restrictions ^{z,r) map S"^ Cq and are quasi-spherical, for 
each {z,r) G Q. As previously, we shall assume $ is 

In order to compute the pullback $*((is|^) using the above techniques, we first rename 
the polar coordinates in Eqn. (^TJ) from {r,6) to (p, C)- Thus we now regard the metric 
parameters U, V, T as functions of the coordinates {z, p, () on the range M x Co of $, and we 
reserve {z, r, 6) for coordinates on the domain Q x S"^. 

The map $ may be described using the Lorentz boost and rotation functions w e 



R e C°°(fi,50(3)) via 



$(2;, r, I 



{z, p, = {z, p{z, r, 9), C{z, r, 9)), 



where as before. 



p{z, r,9) = r f{9; w{z, r)) = r (a/1 + |wP + ^'^w), 
C(^, r, ^) = C(^; u'l^, r), r)) = /-^/^(w + W9). 



(42) 



(43) 
(44) 



To compute the pullback $* ((is|p) we use Eqn. 



5p 



and the definitions Eqs. (|22D, (123|), 



ar 

/ + + b~^wf{bti + sixw) 

f + rw^ieti + 9xsi) + rib + w'^9)9^h 



Defining the 3-vector quantities So,to,7 by 



^0 
So 
to 

7 



WR-^R,W 



W . X w 



we similarly find 



and thus 



'0 I b+l' 

b^^Sow + W~^w^z = b~^{w^z — So X w) 
r{eto + 9x So), 



P,z = rf^z = rf9^to + w^^ 



^*{dp) = (/(I + r9^ti) + w'^(3) dr + {rf9'^to + ^^7) dz + rw'^d9. 
Using the identities Eqs. ([TBI) , (^6|) and the analogous 

pC,z = RAj, 

we also have 

$*(rf() = (,.c/r + C^zdz + Cfid9 = p'^RA{(3dr + ^dz + rd9). 



(45) 



(46) 
(47) 
(48) 
(49) 



(50) 



(51) 



(52) 



(53) 
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We denote the pullbacks of the metric functions U, V, F by a tilde, so for example U = 
^*{U) and U{z,r,e) = U{z, p{z,r,9),Ciz,r,e)). Substituting Eqs. into Eqn. (|l]) 

with coordinates (p, () replacing (r, 6) as already mentioned, gives 

^*{dslp) = -2U{f{l + re'^ti) +w'^l3)dzdr 

- 2U{V + rfe^to + w^^) dz^ - 2rUw^de dz 

+ \RA{r de + (3dr + {^ + A'^R^t) dz)\\ (54) 

Now recall that T = T{z, p, () is angular with respect to the (p, () coordinates, and note 
that the pullback of C^F = simplifies using Eqn. (pOD to 

= r'9^A~'R^t, 

where F = $*(F) = r{z, p{z,r,6),({z,r,6)). This shows that the vector A~^R^r is purely 
angular in the (r, 6) coordinates, and the final term of Eqn. (|5^ becomes 

\rde + /3dr + + A~^R'^f)dz\'^. 

Note also by Eqn. (|T|) that 

A-'R^t = QA~'R^f = QA^AQA-'R^f = QA^R^f. (55) 

Introducing the angular vector 

7 := 7 + A-^R^f - U&w = 7 + e(A^i?^f - Uw) (56) 

and noting that w^A~^ = fvo^, the pullback metric becomes 

^*{dslp) = \rde + (3dr + -fdz\'^ ~2Uf{l + re^ti)drdz 

- 2U{V + rfe^to + lU\Qw\^ - fw^R^T) dz\ (57) 

Comparing this metric with the general NQS metric Eqn. (|I]), which may be written in 
3-vector notation as 

ds%Qg = —2u dz {dr + v dz) + \rd6 + /3 rfr + 7 dz\^, (58) 

we obtain the main transformation result for shear-free metrics. 

Proposition 5 Suppose $ G C°°(i7 x S'^,M x Co) for some domain dM? and $ is quasi- 
spherical with respect to the expanding shear-free NQS metric ds'^^p given by Eqn. f|^i[), with 
null coordinate z. Then $ is described by Eqn. (^^ ) with Lorentz boost and rotation functions 
w e C°°(fi,M3)^ g C'°°(n,50(3)), and the pullback ^*{dslp) is given by Eqn. (\^. 
Defining the pullbacks U = $*(f/), V = ^*(y), F = <I>*(F) and the derived vectors sq, si,to,ti 
in terms ofw,R via Eqs. 1^,1^,1^, and 1^,1^,1^, the NQS parameters u,v,l3,'y 
of the metric $*((is|p) are given explicitly by 

u = Uf{l+re'^ti) (59) 

uv = U{V + rfe'^to + |t/|ewp - fw'^R^t) (60) 

13 = vQti + re X si, (61) 

7 = rQto + rexso + A-^R^t - UQw. (62) 

Moreover, ^ is a diffeomorphism if the vector ti{z,r) defined by Eqn. satisfies 

r|ti|<l, W{z,r)en. (63) 
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Proof: Because ds'gp is non-degenerate (by assumption), $ will be a diffeomorphism iff the 
pullback is also non-degenerate, and this holds exactly when u, the coefficient of dz dr in 
^*{dsgp), is non-zero. But f/ 7^ by assumption, and / = a/1 + \w\'^ + > for all 
w & M^^, 9 E S'^, so the condition reduces to 1 -f- rO'^ti > 0. Clearly this holds for all 9 E S'^ 
if and only if Eqn. ( p3| ) is satisfied. All other statements of the proposition follow from 
previous computations. QED. 



IV. EXAMPLES 
A. Spherically symmetric spacetimes 

The metric form 

dsls = -'^Udz{dr + Vdr) + r'^\d9\^, (64) 

with t/, V functions of (2;, r) only, includes the Schwarzschild metric as the special case 
U = 1, ^ = |(1 — 2M/r), M G M. The geometric mass function m = |(1 — g^'^'r^af^h) for the 
general metric (Q) is given by 

2m{z,r) =r{l-2V/U). (65) 

Again switching from (r, 9) to (p, () coordinates in Eqn. ([6^), Proposition ^ describes ds'gg 
in general Lorentz transformed NQS coordinates, with in particular F = and 

U{z, r, 9) = U{z, p{z, r, 9)), V{z, r, 9) = V{z, p{z, r, 9)), (66) 



where p{z,r,9) = rf{9;w{z,r)) = r{b + 9^w), b = ^/l + \w\'^. In general the angular 
dependence of u and uv will be rather complicated, due to the effects of p-dependence of 
U,V in Eqn. (|6D. 

In the special case of the Schwarzschild metric U = 1, 2V = 1 — 2M/p, and the fields 
/3,7 given by Eqs. (|6T|) , (|62|) with F = are both pure i = 1 spin-1 spherical harmonics, 
and u, V satisfy 

u = {b + 9^w){l + r9^ti) = /(I + r9^ti) (67) 
2M 

= - -r. 7^ + ib + 9'^w) {b - 9^w + 2r9^to). (68) 

r(o + 9^ w) 

If the boost w G is constant and R = I, then we obtain the rather simple NQS metric 
parameters 



/? = 0, 7 = -ew, 



u=b + 9^w, 2v = b-9^w- 



which describe the Schwarzschild spacetime in rigidly boosted coordinates. 

Since in general the functions w{z,r) G M.^, R{z,r) G 5*0(3) are arbitrary, subject only 
to smoothness and the size condition Eqn. (|63D, in order to construct challenging exact 
solutions for numerical relativity benchmarking, we might choose w,R in any reasonable 
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manner, keeping w = and i? = / in regions where we wish the solution to remain exphcitly 
equal to the standard Schwarzschild metric. Note that asymptotic decay conditions may also 
be readily determined: for example the natural conditions 

and similarly for R, give bounded u, v and 7 with /3 — as r ^ 00. 



B. Accelerated Minkowski metric 



By moving the base point of the standard future light cone along a timelike curve in 
Minkowski space, we may construct another class of shear-free NQS metrics. The Minkowski 



metric associated with an accelerated null cone foliation was discussed in |]T0[, using a special 
choice of affine parameter on the null rays to determine the radius function. Let z 1-^ 
{p{z),t{z)) G M^'^ be a future-timelike curve in Minkowski space and denote the tangent 
vector by {p,t), with (— ) indicating d/dz. Two possible normalisations for z are f = 1 and 
f = ^1 + Define ^ : M^,! ^ by 



' X' 




p{z) + rO 


T 




t{z) + r 



(70) 



where (2; = t — r, r, 9) are null-polar coordinates and (X, T) are rectangular coordinates on 
M^'^. Note that \1/ maps the future null cone z = const to the future null cone based at 
(x, t) = {p{z),t{z)). The accelerated Minkowski metric ds\j^^ := dT^ + may be 

written 



-{fdz + drf + \pdz + Odr + rd9\^ 

-2(f - e'^p) dz {dr + i(f + e'^p) dz) + \rde + Qpdz\' 



which is a metric in the shear-free NQS form (^Tj), with coefficient functions 

u{z,e) = T{z)-e'^p{z) 



u 

V 

r 



v{z,e) 
r(z,9) : 



{z) + e'piz)) 



Qp{z) 



(71) 



(72) 
(73) 
(74) 



Note that the timelike condition f > \p\ ensures that U, V are both strictly positive. 

Proposition ^ gives the NQS coefficients for the Lorentz-transformed metric ^*{ds\^.j) 
(with ds\]^j written in {z,p,() as before), and we may simplify as follows. The shear /3 is 
given simply by Eqn. (0) and because QA^R^T = QA^R^p, we find that 



7 = r{eto + ex so) + e(-™ + wr^p). 

Since /( = R{w + WO) = RA'^^O, we have 

^ = f - (^p = f - /-^(w; + We)^R^p 

and 



(75) 
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V = \{r + ep) = \{t + f-\w + WefR'p) 
and r = (/ — CC'^)p. This gives immediately that 

u = {l + re'^h){fT-{w + WefR^p), (76) 
and after some computations, 

2v = ^ ^ ^^^^ ((6 - 9^w)f -{w- WefR^p + 2r9^to) . (77) 

Notice that the spherical harmonic decompositions of u,uv contain terms with i = 0,1,2 
whereas /?, 7 are both pure £ = 1. 

The coordinate system constructed in |]10| corresponds to an NQS metric constructed 
from the choices R = I, w = p{z), with proper time normalisation of {p{z),t{z)). The 



metric Eqn. (12) in |jTO| may be transformed to NQS form with NQS parameters f = b, 
u = 1, (3 = 0, 2v = 1 + 2r6'^to, 7 = rGto + rO x sq, where to = W~^p, Sq = [b + l)~^p x p. 
Thus if the acceleration p is non-zero then both f,7 will be unbounded as r 00. The 
transformation between the NU and NQS coordinates amounts to redefining the angular 
variables (the null cones and quasi-spheres are unchanged), and has the effect of moving the 
conformal isometry P (cf Eqn. (13) of [|I^) to the NQS field 7. 

Alternatively, the choice w = R^p, f = b, with R = jp^ijpp'^ —pp^)R gives a Fermi- Walker 
transported spatial frame. In this case we have *Sq = -^R^p x p, sq = 0, to = W~^R^p, 
and u = 1, 2v = 1 + 2r6^to, /5 = 0, 7 = r0to- I am indebted to Andrew Norton for this 
computation. Note that in this case, the parameters (t, p) and {w,R) may be recovered 
from the metric data to{z) by solving ib = btQ, R = (b + \)~^R{wtl — tQW^) and p = Rw, 
with initial conditions w{zo) = 0, R{zo) = I and p{zo) = corresponding to an initial frame 
at rest. 



C. Robinson- Trautman metrics 



It was shown by Robinson and Trautman |jTl]] that vacuum spacetimes which contain a 
null geodesic congruence which is hypersurf ace-orthogonal, expanding and shear-free have 
particularly simple structure. A coordinate transformation [|I2| brings such metrics to the 
NQS form 

dsjiT = -2Udz {dr + i(Aot/ + U - 2MU-'^ /r) dz) + \rde - U^edA^ (78) 

where M G M is constant, Aq = A5'2 is the standard metric Laplacian on S"^ and U = U{z, 9) 
is independent of r. The vacuum Einstein equations are satisfied by ds\rp if U satisfies the 
nonlinear parabolic equation 

dU 

12M — + U'^AoK = (79) 
oz 

where M 7^ and K = f/^(Ao log f/ + 1) is the Gauss curvature of the metric (/'"^dsQ confor- 
mal to the standard metric ds^ on 5*^. If M = then the metric reduces to the accelerated 
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Minkowski metrics considered above [0. A global existence theorem for Eqn. (^) has been 
given by Chrusciel 

The RT metric has NQS parameters U, V = ^{AqU + U + 2MU-^/r), /? = and 
r = —U^, and after changing Eqn. ([7^) over to {z,p,() coordinates, Proposition || describes 
the effect of a general Lorentz deformation of the coordinates. In particular, since T{z, r, 9) = 
— f/J(2;, T, ^)), using the identity 

Ue = DeUiz, () = r%<;RAQ 
and Eqn. (|55D , we may simplify terms involving F: 

fw'^R^V = w'^A-^R^f = -fw^il^g. 
Thus in the RT case Eqs. (|5I1|),(|5^ become 



uv = U{V + rfe'^to + lU\Qw\^ - fw'^U^g), (80) 
7 = rQto + r9xso- {fU)% (81) 

and u = fU{l + rO'^h), (3 = rQti + r9 x si. 

Note that if the Lorentz deformation preserves j3 = then w = w{z), R = R{z) and 
u = fU, and 7 will remain independent of r only if Sq = to = 0. This requires that w, R are 
constant, and the transformed metric will again be in the explicit RT form of Eqn. (|78D . This 
global Lorentz transformation may be used to normalise to zero the i = 1 spherical harmonic 
components of lim^^oo u{z, 9) (or equivalently, of lim^^oo 7(2, 9)) — this transformation may 
be interpreted as defining an asymptotic rest frame for the RT spacetime fl^, [|l3|] . 



V. DISCUSSION 

In the case of vanishing shear aj^p = (and assuming non-zero expansion pjvp 7^ and 
spherical sections), we have seen that the null hypersurfaces are isometric to the standard 
cone Co (Proposition H), and the residual freedom in the NQS gauge consists precisely of 
a Lorentz transformation at each quasi-sphere. The transformed metric has NQS shear f3 
consisting purely of £ = 1 spherical harmonics, and conversely, if (3 is pure i = 1 then there 
is an inverse quasi-spherical map which transforms the metric into NQS form with (3 = 0. 
Thus the i = 1 spherical harmonic components of the NQS shear (3 are pure gauge. 

Generalised Lorentz transformations preserving the condition (3 = have parameters 
{w, R) depending only on z, since Eqn. (pl| ) combined with Eqs. (|37D,(|38D show that w^r and 
R r must vanish. This remaining gauge freedom may be used to set the i = 1 components 
of 7 to zero at one fixed radius tq as follows. The six i = 1 coefficients of the terms 
Q{A^R^r — Uw) of 7 (cf. Eqn. ([6^)) form a nonlinear functional of w, R, so Lemma ^ may 
be used to solve for w 2, R ^ from so,to, giving a system of ordinary differential equations 

-^{w{z,ro),R{z,ro)) = F{z,ro;w,R), 
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where F{z, tq; w, R) is linearly bounded in w. Consequently there exists a solution which is 
global in z, which in turn ensures (after applying the resulting Lorentz transformation) that 
7(2;, ro) has vanishing i = 1 components at each z. 

Alternatively, it might be possible to use the gauge freedom w{z),R{z) to normalise 
the i = 1 components of u{z, tq) using Eqn. (pQ]), since f = b + is pure i = 0,1 and 
ti = by the condition (3 = 0. Note that this remaining freedom is similar to that available 
in the Bondi and Newman-Unti gauges. In any case, it is a plausible conjecture that the 
gauge freedom remaining in the general NQS metric (|l]) may be used to eliminate the £ = 1 
components of (3, and that the freedom remains to make a rigid Lorentz transformation on 
each null hypersurface. 

The interpretation of £ = 1 spherical harmonic components as gauge terms has also been 
noted in the construction of the Regge-Wheeler-Zerilli equations for linearised perturbations 



of Schwarzschild [T^, 0]. The gauge- invariant quantities satisfying the RWZ equa- 
tions are constructed from i > 2 components of the metric perturbations. Furthermore, 
one quantity constructed from the i = 1 components represents (non-dynamic) angular 



momentum [|T8| , arising from the Kerr perturbation of the Schwarzschild metric — this 
quantity corresponds to the linearised limit of the odd (rotational) £ = 1 component of 
dz{P/r) — dri^j/r), and vanishes for the pure gauge variations constructed above. 
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APPENDIX A: SHEAR-FREE SPACETIMES 

Let A/" be a null hypersurface in some spacetime, with induced (degenerate) metric dsj^ = 
g^f. An adapted null frame on A/" is a pair of vector fields (/,m) where / is a degeneracy 
vector for dsjj-, m G TA/" C, and (/, m + m,i{m — fh)) form a real basis for TM, such that 

gMijn, m) = g^{l, I) = gj^{l, m) = 0, gf/im, m) = 1. 

Using V to denote the ambient spacetime covariant derivative, we define the shear and 
expansion of dsjj- (with respect to the null adapted frame {l,m)) by 

shear = a^p = -giVmrnJ), (Al) 
expansion = p^p = —g(Vm^J)- (A2) 

Although we use the Newman-Penrose notation, the importance of the shear and expansion 
of a null geodesic congruence was known prior to — see |jl5|] for example. 

Lemma 6 Let {l,m) be an adapted null frame for , then the shear and expansion depend 
only on {I, m) and ds'j^. In particular we have 
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(Tnp = -gjvirn, [l,m]) (A3) 
PNP = -HqA^, ^]) + 9Mijn, [/, m])), (A4) 

where is the Lie bracket. 

Proof: The identities ( [A3D , ( [AID are easily verified since [/, m], being the Lie bracket of vector 
fields tangent to the hypersurface A/", is again tangent to M . QED. 

Lemma 7 Suppose (/', m') is a null adapted frame which presents the same orientation of 
M and the null generators as {l,m). There are real functions a,X,fi& C°^(A/') such that 

m' = e'^m + al, I' = e^, (A5) 

and the shear and expansion satisfy 

^T^P = e^+2^V^P (A6) 
Pnp = (i^PNP- (A7) 

Consequently the conditions "pnp 7^ 0, aNP = everywhere onAf" are independent of the 
choice of adapted null frame, and we may consider ctnp/ Pnp as a section of a spin-2 complex 
line bundle over M . 

Proof: The representation (|A5|) for the frame change follows directly from the orientation 
and orthogonality conditions. The formula 

[m\l'] = e'^+*^([m,/] + D^pl - i DiXm) + ef'iaDifi - Dia) I 

where Di,Dm are the directional derivative operators, leads directly to Eqs. (|A6|) ,( [A7|) . 
QED. 

Using the foliation of J\f by null generating curves, we may introduce adapted coordinates 
(p, x^, x'^) by requiring (x^, x^) to be constant along the null generators, and then allowing p 
to be any parameterisation of the null generators. In such coordinates the metric becomes 

djSj^ — ho^iydx dx , 

where the indices a, b have range 3, 4 and hat = hab{p, x^, x^). A natural choice of null frame 
is I = dp and m = m'^da, where dp, 83, ^4 are the coordinate tangent vectors. Introducing the 
cotangent vector ruadx"-, where the rua, a = 3, 4 are defined by the requirements m°'ma = 0, 
fff'TRa = 1, we have 

hah = mafrib + rhamh. 

Direct computation using [I, m] = dp{m"')da gives the following expressions for the shear and 
expansion with respect to the coordinate-based null framing (/,m): 

aNP = m^dp^ma) 

= Im^'m'dphab (A8) 
Pnp = l{m''dp{ma) + m^dpirha)) 

= ^h-'dphab, (A9) 
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where = [/la^]"^ = m°'fh^ + rh°-m!'. 

If M is shear-free and expanding then the metric on M may be brought into exphcitly 
NQS form. It should be possible to extend this result to allow some non-zero shear, but the 
proof will be considerably more difficult. 

Proposition 8 Suppose M is a null 3-manifold with everywhere vanishing shear and non- 
zero expansion, and having spatial cross-sections which are topological spheres. Then there 
exist polar coordinates (r, on A/" ^ M x S*^ such that ds\ = r'^^d'd'^ + sin^ d(^'^) . 

Note that the following argument may be easily adapted in case the spatial sections are not 
spheres. 

Proof: Let I = dp, m = m"'da be a coordinate-based null frame for J\f. By Eqn. ( [X^ ) and 
the shear-free condition we have 

Trfm^dphab = 0. 

Now dphab niay be decomposed 

dphab = A^mafrib + rharrib) + Brriamb + Brhafnb, 

where A is real-valued and B is complex-valued. The shear-free condition shows that B = 
and the resulting equation dphab = Ahab niay be integrated along each null generator to give 

hab{p,x) = exp (^j A{s,x)ds^ h^l^ix), 

where h^^i^) = hab{po{x),x) is a fixed metric on S*^. Now the Riemann Uniformisation 
Theorem [|19| shows there is a diffeomorphism $ : S*^ — S*^, x ^— {'d{x), ^{x)), and a function 
e C°°(S'^]R+) such that h^^dx^dx^ = (l)'^{x)^*{d'd'^ + sin"^ dif"^) . Using the coordinates 
{■}), if) to label the null generators gives the representation 



dslf = exp (^j^ ^^(d^'^ + sin^'i? V)- 



Now define the positive function r = r{p,-^,if) by r = exp Since 4pivp = 

h^^dphab = 2A, it follows that 

dr 

— = rpNp > 0, 

hence r is a valid coordinate, and takes the required form in the coordinates (r, t?, 
QED. 
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